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ABSTRACT 

One puzzle in understanding how stars form in clusters is the source of mass - is all of the mass in 
place before the first stars are born, or is there an extended period when the cluster accretes material 
which can continuously fuel the star formation process? We use a multi-line spectral survey of the 
southern filament associated with the Serpens South embedded cluster-forming region in order to 
determine if mass is accreting from the filament onto the cluster, and whether the accretion rate is 
significant. Our analysis suggests that material is flowing along the filament's long axis at a rate of ~ 
30 Mg Myr -1 (inferred from the N2H + velocity gradient along the filament), and radially contracting 
onto the filament at ~ 130 M© Myr~ 4 (inferred from HNC self-absorption). These accretion rates are 
sufficient to supply mass to the central cluster at a similar rate to the current star formation rate in 
the cluster. Filamentary accretion flows may therefore be very important in the ongoing evolution of 
this cluster. 



1. INTRODUCTION 

Understanding how clusters of stars form is a topic of 
considerable interest, since most st ars appear to form 
within clustered en vironments (e. g.. lLada fc L ada 2003; 
iBressert et all [20Toh . As noted in lMversI (|2011afl . young 
stellar clusters are tightly packed with stars in their 
centres; the initial gas mass density inferred to create 
these stars is significantly higher than observed forma- 
tion conditions. One simple solution to this problem is 
to feed in additional mass as the cluster is forming. Anal- 
ysis of competitive a ccretion simulatio ns indeed show 
that this can occur (jSmith et all I2009T ). Observation- 
ally, young embedded cluster-forming systems are often 
found to be associated w ith multiple fi laments, in a hub- 
spoke like geometry (e.g., Myers 2009a). Recent Herschel 
observations have shown that filaments are ubiquitous 
in star-forming regions, and furthermore that prestellar 
cores and protostars tends to be associated only with 
those filaments w hich are gravitationally supercritical 
(jAndre et all l2010h . Filaments therefore appear to be 
the prime candidate for providing a continuous source of 
material to be accreted onto forming clusters. Few obser- 
vations yet exist, however, to test this idea, as most have 
focussed on t he column density structur e of filaments and 
clusters (e.g.. lArzoumanian et alll2~011l ). 

Several decades ago, low resolution CO observations 
showed that some filaments have velocity gradients 
running along the ir leng th (on ~ 10 pc sc ales, e.g., 
iMcCutcheon et all Il982t iDobashi et al 11992ft . Addi- 
tionally, some filaments had hints of a velocity gra- 
dient across their short axis (e.g., lUchida et al 119901 
IDobashi et all ["l992). Higher resolution observations of 
dense g as species, however, hav e been lacking until re- 
cently. lHacar fcTa falla (20TH) observed the low-mass 
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star-forming region L1517 in Taurus, and showed it con- 
tains several coherent, roughly thermal filaments which 
are well-fit by an isothermal cylinder model. These fil- 
aments have small modulations in their (column) den- 
sity and velocity profiles along the long axis, suggesting 
core accretion from the surrounding filament material. 
Hacar et al (in prep) conducted a similar study of the 
filaments in B213 in Taurus, finding evidence that the 
singular filament visible in column density maps is ac- 
tually composed of many distinct velocity components 
of gas, each of which has a coherent and thermal ve- 
locity profile. In more massive filaments, spectral ob- 
servati ons tend to indica te much more dynamic gas mo- 
tions. iMiettinenl (|2012ft made pointed observations of 
C 18 and C 17 dense clumps embedded in the filamen- 
tary IRDC G304.74+01.32 and find strong signs of in- 
fall and strongly non-thermal linewidths. If the non- 
thermal linewidths can be interpreted as providing addi- 
tional support, then the clumps and fila ment may be in 
yirial equilibrium. In the DR21 filament, Sch neider et all 
(2010) find signatures of infalling motion over both the 
filament and clump scales, and kinematic signatures con- 
sistent with the filament forming from converging flows. 

In order to study the question of filament accretion 
onto clusters, we present kinematic observations of the 
Serpens South embedded cluster and surrounding fila- 
ments. To the best of our knowledge, this is the first 
study to obtain sufficient complementary measurements 
of the gas velocity to provide reasonable certainty in at- 
tributing the motions to filamentary accretion. 

The Serpens South cluster was first discovered in 2008 
by iGutermuth et all from Spitzer data. In this discov- 
ery paper, the cluster is clearly embedded in a dense 
filamen tary structure wh i ch sh ows up in absorption at 
8 (j,m. IGutermuth et all (|2008ft identify approximately 
90 YSOs, and find that the ratio of class II to I sources 
in the region is low, around 0.7, which is much lower than 
the median value of 3.7 in nearby cluster- forming regions 
(|Gutermuth et all |2009 i). and suggests that the Serpens 
South embedded cluster is among the youngest exam- 
ples known. As such, Serpens South may be a particu- 
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larly good target to search for accretion flows in filaments 
and onto the central cluster, since there will have been 
less time for stellar feedback to disrupt or complicate the 
cluster-forming gas mot ions. Based on their Spitzer data, 
iGutermuth et all (|2008| ) estimate that stars are forming 
in the central cluster at a rate of 9 M^JVlyr -1 , although 
more recent work by iMaurv et all (|2011l ) using a combi- 
nation of Spitzer, Herschel, and MAMBO 1.2 mm dust 
continuum observations suggest that the rate may be sev- 
eral times smaller. Much of this discrepancy is due to the 
timescales assumed for protostellar evolution. 

Additional recent observations of Serpens South in- 
clude near-IR polarization measurements suggesting a 
relatively uniform magnetic field a ligned perpendicula r 
to the main filamentary structure (Sugitani et 
and a catalog of protostars associated with outflows sug- 
gesting that the energy injected in the central embed- 
ded cluster by the outflows may be eno ugh to maintain 
the c urrent level of non-thermal motion dNakamura et all 
2011). X-ray observations from Chandra combined with 
the Spitzer data show that the YSOs have similar prop- 
erties to those in Serpens Main and NGC 1333 (e.g, 
iWinston et all [20Tol) . with analysis suggesting that the 
class III sources are likely a by-product of early disk- 
stripping from inter-cluster dynamics, rather than indi- 
cating a population older than the class II sources (Win- 
ston et al, in prep). 

We obtained spectral line maps with MopraB across 
^100 arcmin 2 of the main filament in Serpens South in 
a variety of molecular tracers (discussed in the next sec- 
tion) in order to measure the motion of the dense gas 
in the filaments. We combine this data with measure- 
ments of the column density distribution (Gutermuth et 
al, i n prep) obtained w ith ASTE using the AzTEC cam- 
era (| Wilson et all 12008) to fully interpret the dynamical 
state of the filaments. Using these data, we find strong 
evidence that filamentary accretion flows are playing an 
important role in the ongoing evolution of the Serpens 
South cluster. 

In Section 2, we discuss our Mopra observations, and 
the basic filamentary structure seen in Section 3. In Sec- 
tions 4 and 5, we discuss analyze the spectral data and 
estimate the filamentary accretion rates. In Section 6, 
we perform some simple analysis of the column density 
structure of the filament, then discuss our main results 
in Section 7, and conclude in Section 8. 

2. OBSERVATIONS 
2.1. Data 

Our observations were made with the ATNF Mopra 
22 m telescope in Australia, in September 2008, using the 
single pixel broadband spectrometer, MOPS0. MOPS al- 
lows for multiple spectral windows to be observed simul- 
taneously, increasing the mapping efficiency for multiple 
spectral lines; the spectral resolution in the zoom mode 
is ^0.108 km s _1 for the 90 GHz receiver. We used eight 
zoom windows to focus on different molecules of interest 

5 The Mopra radio telescope is part of the Australia Telescope 
National Facility which is funded by the Commonwealth of Aus- 
tralia for operation as a National Facility managed by CSIRO. 

6 The University of New South Wales Digital Filter Bank used 
for the observations with the Mopra Telescope was provided with 
support from the Australian Research Council 



- the (1-0) transitions of N 2 H+, HNC, HCN, H 13 CN, 
HCO+, H 13 CO+, plus SiO(2-l) and NH 2 D(1,1). Most 
of these transitions trace dense gas; HCO + is sensitive 
to outflow s, while SiO traces pr i marily shocked emission 
(see, e.g.. iSanhueza et all l20ll IBergin fc TafalTal 120071: 
iVasvunina et all I2011L for a discussion on the physical 
conditions traced by these molecules). 

We made on-the-fly maps, scanning along R.A., and 
Dec, with 11" spacing between scan rows, and a scan 
speed of 3.71" s . Square maps of 5.3' per side (with 
the inner 5' fully sampled) were made and later placed 
on a single grid using Gridzilla. An OFF position of 
(18 h 16 m 56 s , +03°02'19") was observed roughly every 
three minutes, and a measurement of the system temper- 
ature taken every thirteen minutes. Every map required 
roughly an hour of integration time for each of the R.A. 
and Dec. scans. We used Livedata 3 to apply the band- 
pass calibration and remove a low order polynomial base- 
line. We then created uniformly gridded spectral cubes 
for each waveband using GridzillsQ, using a final pixel 
size of 15"; Mopr a's spatial resolution is ^40" at 3 mm 
(|Ladd et alll2005h . We converted the spectra from T* to 
T m b units using an antenna efficiency of 0.49 (|Ladd et all 
2005). Due to the relatively large separation between the 
adopted OFF position and our mapping area, with this 
initial reduction, we found some spectra still showed sig- 
nificant ripples in their baseline. In these cases, CLASf-O 
was used to fit and remove a higher order polynomial 
baseline. 

Figure [T] shows the region we mapped with Mopra, 
compared to a dust emission map recently obtained 
with AzTEC/ ASTE by Gutermuth et al (in prep). In 
Figure [U the greyscale image and contours show the 
dust emission, while the blue contour shows the area 
we mapped with Mopra. The yellow diamonds, or- 
ange open circles, and red filled circles show the YSOs 
id entified by Gutermuth in a re-analysis of the data 
in IGutermuth et all (|2008l ) which now includes Spitzer 
24 urn data, and ado pts the classification scheme of 
IGutermuth et al (2009). Clearly, the Mopra observa- 
tions span the main filamentary structure in Serpens 
South. In total, we have more than 1500 spectra at 
each of the 8 wavebands, and cover an area of roughly 
100 arcmin 2 . The median value of the rms noise per 
channel in the spectra are ~0.34 K km s _1 for all eight 
molecular species. Table 1 summarizes the central fre- 
quency adopted for each line and the associated refer- 
ence. For lines with hyperfine structure, the transition 
associated with the adopted frequency is also listed. 

2.2. Line Fitting 

For the emission lines without hyperfine structure 
(HCO+, H 13 CO+, and SiOfl we fit each spectrum 
with a single Gaussian, using the least squ ares min- 
imization code MPFIT dMarkwardtJ 120081 see also 
|http : //purl . com/net/mpf it[ ) to constrain the fit pa- 

7 Livedata and Gridzilla are ATNF software for on-the- 
fly maps, written in AIPSH — h They are available at 
www . atnf . csiro . au/ computing/ soft ware /livedata/ index .html. 

8 CLASS is available at www.iram.fr/IRAMFR/GILDAS. 

9 H 13 CO+ has b een s hown to have hyperfine structure by 
ISchmid-Burgk et a 3 l l200l) . but the lines are spaced too closely 
to be discernable in our data, and so we treat the line as having 
no hyperfine structure. 
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rameters to a reasonable range and calculate the errors 
associated with each fitting parameter. We similarly fit 
NH2D with a single Gaussian, as the line was sufficiently 
weak that none of the hyperfine satellites were detected 
(in most cases, not even the central component was visi- 
ble). For these four transitions, the parameters fit are the 
velocity centroid, velocity dispersion (in Gaussian sigma 
units), and line peak. 

For the emission lines with hyperfine structure (N2H" 1 ", 
HNC, HCN, and H 13 CN), we substitute the Gaussian 
model with a single uniform slab model with the ap- 
propriate hyperfine em ission information (adapted from 
iRosolowskv et all I2008Q . and again use MPFIT to find 
the best fitting line parameters - in this case, the ve- 
locity centroid, velocity dispersion (in Gaussian sigma 
units), optical depth, and excitation temperature. 

For all of the results discussed in this paper, we only 
analyze spectra which are above a signal to noise level 
of 3.5 - i.e., we require an integrated intensity to be 
3.5 times that expected from the noise level (the spec- 
tral rms times the square root of the number of velocity 
channels), and also a minimum of two spectral channels 
showing emission about 3.5 times the spectral rms (us- 
ing two channels instead of one decreases the likelihood 
of the emission being a noise spike). Calculation of the 
integrated intensity is discussed in Section 3. 

Over nearly the entire area mapped, most of the 
emission lines detected show strong indications of self- 
absorption; only N2H + , H 13 CO + , and, where detected, 
NH2D appears single-peaked, while SiO has no un- 
ambiguous detections. Even the N2H + lines may be 
mildly self-absorbed in some locations, but our signal 
to noise and spectral resolution are insufficient to de- 
termine this definitively. Figure [5] shows the spec- 
trum of all eight emission lines near the cluster centre, 
at the location where the N2H+ emission is strongest, 
(-18' l 30 m 01 s , -2°03'12"). In all cases, the spectra have 
been summed over the nearest four by four pixels (60") 
to increase the signal to noise. The red vertical lines 
indicate the location of the lines based on the centroid 
velocity fit to the summed N2H + spectrum. Hyperfine 
components are indicated, along with their relative inten- 
sities expected for the optically thin case. As the figure 
clearly illustrates, many of the emission lines have high 
levels of self- absorption. We therefore focus primarily on 
the line fit parameters for the N2H 4 " and H 13 CO + emis- 
sion. Note that since there are seven N2H + hyperfine 
components with different optical depths, the N2H+ line 
widths and centroids fit tend to be satisfactory, even if a 
small amount of self-absorption is present. 

3. BASIC FILAMENT STRUCTURE: INTEGRATED 
INTENSITY 

We first calculate the integrated intensity for the five 
molecular lines that do not exhibit self-absorption and 
hence are more likely to be optically thin to show the col- 
umn density structure of the dense gas. The full MOPS 
bandwidths are very large (~200 km s _1 ); in order to 
minimize the contribution of noise to the calculation, 
we restrict our integration to velocity channels close to 
where emission is expected. The velocity range is suffi- 
cient to extend slightly beyond the maximum observed 
line extent; the range adopted for each emission line is 
given in Table 1. The velocity centre for integration is, 



in order of preference, as the signal to noise level allows: 
the centroid velocity fit to the particular line, the cen- 
troid velocity fit to the N2H+ lineF 7 !. or the approximate 
cloud-wide centroid velocity, 7.5 km s . Figure |3] shows 
the integrated intensity maps for N2H + , H 13 CO + , and 
NH2D, i.e., the brightest species without self-absorption, 
along with SiO, which lacks any clear detections. The 
typical uncertainty is given in Table 1. 

Figure [3] shows that both the N 2 H+ and H 13 CO+ inte- 
grated intensities have a similar structure to that seen 
in dust emission. The integrated intensity for NH2D 
is much fainter, but appears to follow a similar trend. 
All three of these molecular species are expected to be 
good tracers of dense gas (see, e.g., the di scussion in 
iSanhueza et~al I20H iBergin fc Tafahal [20071 ) . which is 
confirmed by their good correspondence to the dust emis- 
sion. H 13 CN (not shown in Figure El) is also expected t o 
be a good dense gas tracer (e.g., Hily-Blant ct al 2010), 
however, our data has too poor a signal to noise level to 
see this definitively. Summing our data over many pix- 
els shows hints that the central filament is brighter than 
the surrounding gas, although even this is only seen at a 
marginal level. The H 13 CN spectrum shown in Figure [2] 
is one of the brightest spectra resulting from summing 
the data over four by four pixels squared. 

SiO is not detected with certainty any where in the 
map. As a tracer of shocked gas (see e.g., iCaselli et all 
119971 and references therein), its emission structure could 
be expected to differ from the other four species. The 
general lack of emission may indicate a lack of strong 
(highly supersonic) shocks, however, with the noise lev- 
els in our observations, we cannot rule out weaker shocks. 

3.1. Filament Geometry 

In Figures [T] and El differing filament shapes or ge- 
ometries are visible north and south of the main cluster. 
The southern filament is dominated by a single nearly 
straight filament extending to the southeast, with sev- 
eral fainter filaments (e.g., pointing due south) appar- 
ent primarily in the dust emission map. In the north, 
however, the filamentary structure is more complex. A 
careful visual examination of both the dust emission 
and molecular data to the northwest of the main clus- 
ter shows a seemingly complex structure suggestive of 
overlap along the projected line of sight. Intertwined fil- 
aments have recently been identified through a careful 
analysis of the quiescent B213 filament in Taurus (Hacar 
et al i n prep). Other observations of filaments, s uch as 
Orion (j Johnstone fc Ballvl [19991 iBallv et alll987t) show 
morphologies suggestive of similar phenomena in higher- 
mass systems. Unlike B213, the large linewidths in the 
northern filament of Serpens South make the separation 
of velocity components difficult if not impossible, and for 
this paper, we restrict our analysis to the simpler south- 
ern filament. 

We define the central or peak ridge of the southern fila- 
ment by searching for the location of peak emission along 
every horizontal cut through the filament. For this de- 
tection, we use the dust continuum map; while the N 2 H+ 
integrated intensity map shows a similar distribution, the 
dust map has a higher signal to noise level. This filament 

10 N.B.: the N2H + line is the brightest and most frequently 
detected of the five species without self-absorption. 
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definition is shown in Figure 2] For comparison, we also 
show a similar outline for the more complex northern fil- 
aments. For our analysis of the southern filament, we 
consider the data within V of the peak ridge (i.e., three 
15" pixels on either side of the peak, or a diameter of 
0.16 pc), and along a 0.33 pc extent. This area covers 
all of the strongest emission present in the southern fil- 
ament and excludes any gas potentially associated with 
the much fainter filament to the west. Within this area, 
the southern filament has a mass of ~20 M© , based on 
the dust emission ma r£J IMaurv et all lj20"H estimate a 
mass of 610 M© over a 15' by 25' area (i.e., several times 
larger than the entire extent of our Mopra observations) , 
based on 1.2 mm MAMBO dust continuum dat a. For 
both our estimate and that in IMaurv et~al (|201lD , a dis- 
tance of 260 pc (e.g., iStraizvs e t al 2003J), is assumed, 
based on extinction estimates in the Aquila rift. 

There is, however, some debate over the true distance 
to Serpens South. Recent VLBI measurements of a YSO 
in the Serpens main core sugge st that the dist ance to 
Serpens main is roug hly 415 pc (IDzib et aJll20loh . Com- 
plicating the picture is the fact that the gas velocities 
appear similar for Serpens main, Serpens So uth, and 
Aquila; see the discussion in IMaurv et all (poTTh for more 
details. Since most existing analyses of Serpens South as- 
sume the closer distance o f 260 pc Ce.g..lGutermuth et all 
200$ IMaurv et all 120111 : lArzoumanian et all 120111 ). we 
also adopt that value here, but discuss in Section 7 how 
our results would differ if we instead assumed 415 pc. 

4. GLOBAL VELOCITY GRADIENT 

Figure H] shows the centroid velocity fit to our N 2 H+ 
spectra (where the signal to noise ratio is 3.5 or more). A 
clear velocity gradient is visible along the southern fila- 
ment (shown by the red dashed line) , while the northern 
filaments have a structured, but more complex velocity 
field. If the southern filament is oriented with its south- 
east part further from us and the cluster centre closer 
to us, then an accelerating accretion flow would produce 
a velocity gradient similar to the one observed. The in- 
fall line profiles discussed in the next section lend further 
weight to this choice of filament orientation. 

We measure the velocity gradient for t he southern fila- 
ment following the procedure outlined in lGoodman et all 
(fl99l . We assume the N2H 4 " centroid velocities follow 
a simple linear form: 



AAc 



BA8 = C 



(1) 



where Aa and AS are the changes in right ascension 
and declination across the map. We found the best- 
fitting values to t he constants A and B using MPFIT 
(Markwardt 2008), and then converted all of the angu- 
lar velocity gradients to physical scale assuming a dis- 
tance of 260 pc. This yields an overall velocity gradient 
of 1.4 ± 0.2 km s _1 pc _1 , at an angle of ~ 20° west of 
north. Note that while the H 13 CO + centroid velocities 
(not shown) agree quite well with the N2H + centroid ve- 
locities, fewer of the spectra have a signal to noise level 

We assume a temperature of 1 5 K and a dust opacity of 
0.0114 cm g~ - 1 (e.g. Enoch ct al 2007); recent ammonia measure- 
ments suggest a dense gas temperature around 10-15 K in Serpens 
South, with the southern filament at ~11-12 K further from the 
central cluster and ~12-14K approaching the cluster (Friesen et al, 
in prep). 



above 3.5, which leads to a less reliable measurement of 
the velocity gradient. 

We can estimate the accretion rate implied by this 
velocity gradient by using a simple cylindrical model, 
shown in Figure[SJ The cylinder has a mass M, length L, 
radius r, an inclination to the plane of the sky of angle 
a, and true motions of velocities Vj| along the filament 
long axis, and V±_ inf ailing into the short axis (this final 
term will be of interest in the following section). The 

accretion rate, M\\, along the filament onto the central 
cluster is therefore the velocity along the filament times 
the density and the area perpendicular to the flow, or 



M|| = V\\ 



M 



(irr 2 ) 
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which simplifies to M11 = Vji . Due to projection ef- 
fects, we observe 

L obs — Lcos(a) and Vj|, bs = V\\ sin(a) (3) 

with Vn b s — VVii (, s L f, s . Therefore, the accretion rate 
is given by 
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For the southern filament, we estimate a mass of 
20.3 M along a length of 0.33 pc and width of 0.08 pc 
(see Section 3.1), which corresponds to an accretion rate 
of 28 Mq Myr -1 for tan(a) = 1. Note that since we ob- 
serve both a gradient along the filament and infall onto 
the filament (discussed in the next section), this already 
constrains a not to lie close to either of the possible ex- 
treme values (0° for lying parallel to the plane of the sky 
and 90° for lying directly along the line of sight). We will 
furthermore show that the accretion rates estimated from 
both the motion along and across the filament are rela- 
tively high, which suggests that a is likely close enough 
to 45°, so that the correction factor for each measure is 
only several times unity. 

The velocity gradient used for the accretion rate es- 
timate is also similar to values measured i n oth er star- 
forming regions - for example, Ball v et all (|1987l ) find a 
gradient of ~0.7 km s _1 pc -1 along several pc of the 
Orion integral shaped filament; I Schneider et all (j2010f ) 
find gradients of ~0. 8-2. 3 km s" 1 pc" 1 along the DR21 
filament, and [Kirk et al (|2010D find a general relation- 
ship in non-filamentary molecular gas structures be- 
tween size and velocity gradient suggesting values around 
1 km s _1 pc -1 for size scales of a few tenths of a parsec. 
The similarity of the measured gradient further suggests 
that projection effects are not extreme for the filament, 
if Serpens South is not atypical. 

If the gradient were instead interpreted as rotation, 
then the total rotational energy would be 



^-I(^-)(W)" 



(5) 



where the second term is the moment of inertia for a 
cylinder rotating longways around its end, and the final 
term represents the angular velocity squared. For the 
southern filament, this yields E rot ~ 0.7 M 
The gravitational binding energy is given by 



km 2 s 2 . 
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where M c i ust is the mass of the central cluster, roughly 
100 M Q , calculated in the same manner as the filament 
mass (see Section 3.1). Thus, the gravitational binding 
energy is a factor of more than ten higher, at E grav ~ 
26 M Q km 2 s -2 . With a freefall time of just ~0.3 Myr, 
several times less than the estimated age of the region, 
even if the filament was initially rotating, infall motions 
would rapidly dominate. 

Finally, the velocity gradient could also be interpreted 
as outflow motion, if the orientation was the reverse of 
that assumed above. In this instance, stellar winds and 
outflows would be the most likely candidate source for 
driving the motions. Using a simple calculation to es- 
timate the momentum and energy associated with such 
outflowing motions yields va lues of 9 M ff km s" 1 and 
2 M km 2 s -2 respectively. INak amura ct al (20ll|) es- 
timate that the total momentum and energy in out- 
flows (measured using CO(3-2) data) in Serpens South is 
~ 8 M© km s" 1 and ~ 65 M© km 2 s _1 respectively, with 
a total mass of outflowing material of only 0.6 M© . While 
in principle, the YSO population could provide sufficient 
energy and momentum to dri ve the observed l arge-s cale 
flow, the outflows observed by INakamura et al (|2011h are 
multi-directional, not highly collimated, as would be re- 
quired. 

5. INFALL ALONG THE LINE OF SIGHT 

We next utilize information from an optically thick 
line showing self-absorption to measure gas accretion 
across (i.e., building up additional mass onto) the fila- 
ment. In the classic case of two converging layers, the 
gas is assumed to have an increasing excitation tem- 
perature toward the centre (where the density is also 
higher); self- absorption in optically thick lines occurs for 
the emission originating from the furthest pieces of gas at 
any given velocity. Therefore, emission from the (blue- 
shifted) layer moving toward the observer is emitted pri- 
marily from the higher excitation, denser part of the slab, 
while the reverse is true for the (red-shifted) laye r mov- 
ing away from the observer (e.g. IMvers et al 1996). Infall 
profiles therefore show brighter peaks blue-ward of the 
peak emission of an optically thin line, with a lower red 
shoulder, and sometimes a central dip. Using simplified 
assumptions about the gas properties, the emission pro- 
file of a line showing such a red-blue asymmetry can be 
used to estimate a characteristic inf all velocity of the gas 
(e.g., lEvanslll999t IMvers et alll2000l ). 

In Serpens South, we observe self-absorption in sev- 
eral of our emission lines: HCO + , HCN, and HNC. Of 
these lines, the HNC is the least self-absorbed, and there- 
fore is the brightest and has the highest signal to noise 
level. A high signal to noise level is r equired to obtain 
reason able model fits to infall spectra - IdeVries fc Mversl 
(2005) find signal to noise levels > 30 are required for 
good results. In order to further maximize our signal to 
noise level, we therefore spatially average spectra along 
the southern filament. As discussed in the previous sec- 
tion, there is a global velocity gradient running along 
the southern filament, which needs to be accounted for 
prior to averaging. We use the N2H+ centroid veloci- 
ties fit at each pixel to first shift all of the HNC spectra 
to a common velocity centre prior to averaging, again 
considering spectra within 1' (or 3 pixels on either side) 
of the filament's peak ridge identified. Figure [6] shows 



the result of this summation for HNC, as well as for the 
two optically thin lines, N2H+ and H 13 CO + . In the case 
of N2H + , only the isolated component is shown for sim- 
plicity. The HNC spectrum shows a brighter blue and 
fainter red peak with strong central dip around the peak 
emission of N2H 4 " and H 13 CO + , in agreement with ex- 
pectations for infall. The HNC (1-0) transition does have 
hyperfine structure, however, this is on a much smaller 
scale than the infall profile seen, as shown in the right 
hand panel of Figure El 

We estimate the infall velocity usin g the Hill5 model 
and line fitting code) described in IdeVries fc Mversl 
2005) . The Hill5 model considers gas with a linearly in- 
creasing excitation temperature toward the centre, with 
each half of the gas movi ng toward the other at a fixed 
velocity. In their analysis. IdeVries fc Mversl (|2005| ) found 
that the Hill5 model tended to provide the best fit to a 
range of both simulated and real spectra with infall mo- 
tions; the inclusion of an increase in the excitation tem- 
perature allowed f or much better ma tch than the classic 
two-slab model of IMvers et all (jl996). The variables in 
the Hill5 model are the infall velocity, Vi n , the (Gaus- 
sian sigma) velocity dispersion, 8V, total optical depth, 
t, system velocity Vi sr , and peak excitation temperature, 
Tpeak- Note that the system velocity is expected to be 
in our case, since the individual HNC spectra were all 
shifted by the N2H + centroid velocity fit at the same 
pixel prior to the summation. Using the Hill5 model, we 
find the best-fitting parameters listed in Table 2 under 
'Fit-1'. 

In order to estimate the fitting errors, we ran the fit- 
ter on 1000 realizations of the spectrum with random 
noise added at the same level as observed 12 ! This test 
yielded several unexpected results. In a small number of 
cases (26 of 1000), the 'best' fit clearly got 'stuck' in a 
local minimum of \ 2 space, with unphysical parameters 
(optical depths in excess of 20 and peak excitation tem- 
peratures often around several hundred). These models 
provided visually poor fits to the spectra, and so we dis- 
carded these cases. These fits likely could have been fixed 
by modifying the fitting tolerances in the Hill5 fitter, but 
this seemed unnecessary, given the small overall number. 
Of the remaining 974 fits to the data plus random noise, 
686 were found to have values clustered very close to the 
original fit (these were used to calculate the formal er- 
rors in 'Fit-1'), while the remaining 288 were clustered 
around a different set of fit parameters, listed under 'Fit- 
2' in Table 2. Visually, this other set of fits also are a 
reasonable match to the data, and the distribution of 
X 2 values are similar for both sets of fits. Therefore, 
these second fit parameters appear to be a viable second 
solution to the best-fitting Hill5 model, using only the 
information available from the HNC spectrum. 

We use the additional information available from the 
optically thin spectra (H 13 CO + and possibly N 2 H+) to 
argue that Fit-1 best matches the full suite of data. 
We calculate the corresponding optically thin spectra for 
both Fit-1 and Fit- 2, using the fitted parameters but a 
lower optical depth. For Fit-2, the optically thin spec- 
trum predicted has two distinct velocity components sep- 
arated by more than the linewidth. Effectively, Fit-2 

12 Note that this increases the total noise in the spectrum by a 
factor of V2- 
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represents a model of two independent layers of mate- 
rial moving toward each other. As shown in Figure [6l 
two distinct velocity components are not observed in ei- 
ther of the brightest optically thin tracers, i.e., the N2H + 
isolated hyperfine component or H 13 CO + . This implies 
that Fit-1 is the most consistent with the observations. 

We ran several simple tests with t he online non-LTE 
radia tive transfer calculator, Rade>0 (|Van der Tak et all 
2007) to check that the fitted values are reasonable. We 
assume a n HNC abundance of 6 x 10~ 10 typical of dense 
core gas (|Hilv-Blant et allfeOlOl ). and adopt gas proper- 
ties (density, column density, temperature, total velocity 
dispersion) based on the observations discussed earlier. 
The predicted optical depth ranges from roughly 6 to 
8, for temperatures (and the associated derived densi- 
ties and column densities) between 12 and 15 K. Given 
that a cold dense core HNC abundance is likely an over- 
estimate for a filament, and the uncertainties in various 
parameters used to estimate the density and column den- 
sity from the dust maps, the Hill5 model fit of an optical 
depth of 2.7 for 'Fit-1' is in reasonable agreement; the 
optical depth of 1.6 from 'Fit-2' is less consistent with 
the observational expectations. 

Using the same toy model as shown in Figure [5l the 
accretion rate onto the filament can be estimated from 
the above infall velocity, i.e., the velocity onto the fila- 
ment times the density and filament surface area of the 
short axes, or 



A/_l = V±_ x 

which simplifies to M± = 
effects, we observe 



M 
nr 2 L 



(27rrL) 



(7) 



2Vj_ — . Due to projection 



V in = V±cos(a) 



(8) 



while the observed filament radius / short axis should be 
invariant under rotation. Therefore, the accretion rate 
onto the filament is given by 



Mi = 



rcos(a) 



(9) 



This corresponds to an accretion rate of 130 M© Myr -1 
for cos (a) = 1 and is therefore a lower limit to the true 
rate. Additionally, the infall velocity fit with the HillS 
model is expected to be a lower limit to the true char- 
acteristic infall speed: the infall was modelled as two 
infalling parallel slabs, while in a cylindrical geometry, 
some of the inward motions will appear smaller along 
the line of sight. Finally, note that using the alternative 
Vi n of 0.54 km s _1 from Fit-2 corresponds instead to a 
lower limit of ~ 280 M Myr -1 . 

Both the accretion rate onto the filament and from the 
filament onto the cluster, without accounting for projec- 
tion effects, are within a factor of a fe w of the star for- 
matio n rate in the cluster estimated by iGutermuth et all 
(120091). and s omew hat larger than the rate estimated by 
iMaury et all (|2011| ). Any correction for projection effects 
made to the accretion onto the filament will serve to fur- 
ther increase that estimate, i.e., beyond 130 M Myr -1 , 
while the accretion rate from the filament onto the 



13 Radex is available at http://home.strw.leidenuniv.nl/- 
^moldat a/r adex . html . 



cluster could be larger or smaller than our estimated 
28 M Q Myr -1 after correcting for projection. 

If the filament were projected at an angle of 30 to 60 
degrees relative to the plane of the sky, then the accre- 
tion onto the central cluster measured would be within a 
factor of 0.6 to 1.7 of its true value, while the accretion 
onto the filament measured would be a factor of 1.15 to 
2 of its true value. Based on the combination of filament 
appearance, velocities measured, and accretion rates de- 
rived, we believe that 30° < a < 60° is a reasonable 
range of filament angles. 

5.1. Other Lines 

In addition to the self-absorbed HNC emission, both 
HCO + and all three hyperfine components of HCN also 
show signs of self-absorption. We performed the same 
spectral averaging across the southern filament for these 
emission lines as described for HNC. We additionally 
shift each HCN (and associated H 13 CN) hyperfine com- 
ponent to the velocity it would have if it were the sole 
emission line. Figure [7] shows the resulting spectral line 
profiles. As can be easily seen from the figure, all four of 
the optically thick lines show such strong self- absorption 
that relatively little emission is seen. We attempted to 
model the emission lines with a Hill5 (or other) profile, 
but due to the low signal to noise, even with averaging 
over the filament, we were unsuccessful. 

The degree of self-absorption in HCO + and HCN is 
surprisingly large - the HCO + emission peaks at only 
half the value of the H 13 CO + peak, whereas in an op- 
tically thin environment, the ratio would be somewhere 
around a factor of 40 (for non PDR molecular clouds at 
a galactocentric s eparation, similar to the Sun, of 8 kpc; 
iSavage et all2002ft . In H 13 CN, the line ratios for optically 
thin gas would be 1:5:3 (for the top right, bottom left, 
and bottom right panels in Figure [7|), whereas HCN is 
clearly brightest in the (F=0-1) or first / leftmost hyper- 
fine component. This would be consistent with a picture 
wherein the highest levels of self- absorption occur in the 
optically thickest components, and the least in the opti- 
cally thinnest. 

Recent work suggests that the rate coefficients for HCN 
and HNC differ significantly (D umouchel et aj|2010f) , im- 
plying sensitivity to different density ranges of material, 
and thus explaining why the two do not show similar lev- 
els of self-absorption. Alternatively, studies have shown 
that hyperfine anomalies in HCN(l-O) are relatively com- 
mon in higher column density environments, and tend to 
lead to precisely the tre nd in peak brightnesses observed 
(ILouehnane et all l20l"l . The H 13 CN emission is ex- 
tremely faint, but appears to be roughly consistent with 
optically thin emission with the predicted line strengths 
(with no significant emission detected in the first / left- 
most / F=0-1 component), thi s is also consistent wit h 
the hyperfine anomaly study of ILouehnane et al (|2012f) . 
who find that H 13 CN does not have anomalous ratios in 
sources where HCN does. 

We ran several additional checks to ensure that these 
severely self-absorbed spectra are real, and were not af- 
fected by any systematic effects. For the velocity cen- 
troids, all eight molecules were observed simultaneously 
with MOPS, suggesting that any error with either the 
measured line velocities, or the adopted central frequen- 
cies of each molecule, would have to be present in all 
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of our observations. The Serpens South data presented 
here are part of a larger project which includes four 
other embedded cluster-forming regions. Examining our 
data in the other four regions shows very good agree- 
ment between the different molecular line centroid ve- 
locities where there is optically thin emission, suggest- 
ing that there are no systematic effects with our velocity 
scale. Similarly, the peak line emission in these other 
regions, and the relative values of peak intensity, agree 
well with expectations, suggesting no systematic effects 
there either. Finally, we note that similarly severely self- 
absorbed emission line profiles showing a complete ab- 
sence of a red shoulder have also been observed in Ser- 
pens South with a different telescope and molecular line 
transition. CO (3-2) emission observed with the ASTE 
10 m telescope show a similarly shaped profile in some 
areas of the cluster (see, e.g., the bottom panels in Fig- 
ure 6 of iNakamur a et a mom . 

6. COLUMN DENSITY PROFILE 

6.1. Mass per Unit Length 

In this section, we examine some basic properties of 
the column density structure of the southern filament. 
The simplest property which can be measured for fila- 
ments is their mass per unit length. Assuming that a 
filament can be approximated as an infinitely extended, 
self-gravitating, isothermal cylinder with no magnetic 
support, then the maximum mass per unit length which 
can be in equilibrium is 

2c 2 JG (10) 

M HnetCrit = 16-7(— ^-r) M Q pc" 1 (11) 

(e.g.. I0strik"erlll964t llnutsuka & Mivamal 119971 ). Obser- 
vations of the ammonia (1,1) and (2,2) inversion tran- 
sitions (Friesen et al, in prep) indicate that the kinetic 
temperature across Serpens South is in the range of 10 to 
15 K, implying Mu ne ^ cr i t is between 17 and 25 M Q pc _1 . 
If the observed non-thermal motions were also con- 
sidered to contribute to supporting t he filament (e.g., 
iFiege fc Pudritdl2000atlMiettinenll2012fl , then with a typ- 
ical velocity dispersion of ~ 0.36 km s^ 1 (Gaussian sigma 
units) for N2H + (mean or median value over the southern 
filament), then Mu ne:Cr i t is roughly 80 M Q pc _1 . 

In the southern filament, we find a mass of 20.3 M© 
along 0.33 pc, implying Mu ne — 62 M Q pc _1 , i.e., sev- 
eral times the critical value, unless all of the non-thermal 
motion is contributing thermal-like support. The total 
amount of support available is therefore likely insuffi- 
cient to keep the filament in equilibrium, implying that 
it should be radially contracting. Supercriticality is con- 
sistent with our velocity measuremen ts showing that gas 
is glob ally infalling in the filament. lArzoumanian et all 
(| 2 1 If) similarly find that filaments with higher column 
densities, such as in Aquila, tend to have supercriti- 
cal values of Mu ne . Filaments associated with mas- 
sive star formation tend to be highly supercritical with 
Miine. often ten times la rger than that in nearby regions 
(Henncm ann et all2012| ) . 

6.2. Radial Column Density Profile 



In this section, we analyze the radial column density 
profile of the filament using the millimetre continuum 
emission. We use the peak ridge of the filament dis- 
cussed in Section 3, and for each pixel in the original 
AzTEC / ASTE data, we measure the column density pro- 
file along lines of constant declination. This provides for 
more regular and better sampling than to measure points 
perpendicular to the local filament orientation; we later 
correct our radial separation measurements to account 
for an average filament angle of 45°. 

All ground-based (sub)millimetre measurements are 
insensitive to the largest scale structures, due to the 
observing methods which must be employed. For the 
AzTEC/ ASTE map we analyze, Gutermuth et al (in 
prep) have devised a novel technique which is able to 
recover flux out to about 5', a much larger scale than is 
typical. In summary, the technique iteratively subtracts 
a model of the astronomical signal from the bolometer 
timcstream data and then recomputes the atmospheric 
filtering and regenerates the signal model map. The ini- 
tial model is just the base pipeline-produced map masked 
to a signal to noise of 2, and convergence of the process 
is achieved when no new signal is found in an iteration 
down to a signal to noise of 3. Despite this marked im- 
provement in map reconstruction, the very largest-scale 
structures with shallow or flat radial column density pro- 
files will not be included in the map. Some amount of 
large-scale structure is expected to be present, as Serpens 
South is coincident with a roughly degree wide overden- 
sity of dust extending to the well-known W40 active star- 
fo rming region. The 2MASS -based dust extinction map 
of iRowles fc Froebrichl (2009)03 show that visual extinc- 
tions in excess of 5 to 10 mag are common across this 
area. Some, but not all, of this extinction undoubtedly 
represents diffuse material unassociated with the indi- 
vidual star-forming regions like Serpens South. Any ma- 
terial that is associated with Serpens South, however, 
would effectively result in a missing linear offset to the 
derived column density profile; the Serpens South fila- 
ment is sufficiently compact that the shape of the col- 
umn density profile should be correct. The function we 
fit to the column density profile is dependent on abso- 
lute values, not just the shape, so we limit our analy- 
sis to column densities where the fractional errors due 
to any missed large-scale structure are reasonable. If a 
constant background of Ay = 1 mag is missing from 
the AzTEC/ASTE map, there would be a <30% un- 
derestimate of the inferred column density at fluxes of 
0.067 Jy bm _1 ; we adopt this as our lower limit. An 
additional complication to the analysis arises from a sec- 
ond, fainter filament to the west of the one we analyze. 
In order to prevent this structure from confusing the ra- 
dial column density profile, we additionally exclude pix- 
els which lie westward of the minimum value between the 
two filaments. 

After measuring the radial column density profile along 
the entire peak ridge, we calculate the mean and stan- 
dard deviation at each radial separation, similar to 
lArzoumanian et "al (|2011l ). The mean profile and stan- 
dard deviation of the column density profile of the south- 
ern filament is shown in Figure [5] For comparison, the 

14 These maps are available at 

http : //astro .kent . ac.uk/extinction 
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dotted red line shows the ASTE beam; clearly the fila- 
ment is resolved. 

We compare the filamen t column density profile to the 
same suite of models as in I Arzoumanian et all (|2011h : 



£_(r)=A, 

P [l + (r/R flat )^ 



(12) 



where S p is the column density at radial separation r, 
p c is the central (mass) density, Rfiat characterizes the 
inner flat part of the profile, and 



A n — 



1 



du 



P CO«y_ 00 (1+U 2 )P/ 2 



(13) 



is a constant, and the inclination angle along the line of 
sight, i, is assumed to be for simplicity. 

Assuming that the observed column density profile is 
additionally convolved with the ASTE beamsize, we use 
the MPFIT minimization code ()Markwardtll2008[) to de- 
termine the best fitting parameters. Allowing all three 
free parameters to vary, we find p c — 3.6±2.2 x 10 5 cm~ 3 , 
R flat = 0.010 ± 0.006 pc, and p = 1.9 ± 0.2. Our 
values for Rfi„.t a n d p li e within the range found by 
I Arzoumanian et "al ([20111 ) for all filaments within the 
Aquila Rift (note their p c values were not published), 
including a power law index which is shallower than for 
an isothermal cy linder. It is impor t ant t o note, how- 
ever, that unlike lArzoumanian et a] ()201 ll ) . we are not 
able to make a strong a statement on the power law in- 
dex of the filament. The restrictions discussed above on 
where we are able to measure the radial column density 
profile limi t our analysis to re latively close to the cen- 
tral ridge. lArzoumanian et af s comparable radial col- 
umn density profile extends a factor of nearly 30 larger 
in radius than ours. With the relatively small range 
in radial separations, our ability to constrain the power 
law index is weak. Figure [5] shows also the best fit ob- 
tained using a power law index of 4 (corresponding to 
an isothermal cylinder), demonstrating that an isother- 
mal model also can provide a reasonable fit to our data. 
The temperature required for the isothermal fit is 43 K, 
which is similar to the temperature implied by the non- 
thermal velocity dispersion in N2H+ (around 50 K for 
a = 0.36 km s -1 ). This is consistent with Section 6.1 
where we showed that thermal pressure was insufficient 
to balance t he gravity. 

Recently, Uuvela et all (|2012l) created synthetic obser- 
vations of filaments formed in turbulent simulations, to 
investigate the effect of instrumental noise and changing 
dust properties on the inferred filament column density 
profile and mass per unit length, especially as relates to 
Herschel data. They find that the mass per unit length 
is reasonably robust to these uncertainties, while the pa- 
rameters fit to the filament column density profile can be 
more uncertain in cases of higher noise / poorer resolu- 
tion or a dust opacity which varies with density. 

7. DISCUSSION 

The southern filament in Serpens South has similar 
column density properties to star-forming filaments re - 
cently studied with Herschel ()Arzoumanian et "al [20Tl . 
with a mass per unit length several times larger than the 
critical value for thermal support, and a radial column 
density profile which is consistent with p ~ 1.5 — 2.5. 



In fact, the Serpens Sou th filament is included in the 
lArzoumanian et all (|2011h sample of 'Aquila' filaments, 
although the individual filament locations / properties 
are not given in their paper to allow a more detailed 
comparison. The relatively large velocity dispersions 
seen in dense gas tracers in the southern filament could 
be interpreted as evidence for providing sufficient stabil- 
ity against gravitational collapse, however, our more de- 
tailed kinematic analysis suggests that the filament is in 
fact not in equilibrium. We find evidence for mass flows 
both along the filament onto the central cluster and onto 
the filament from the less dense surroundings, with ac- 
cretion rates of order 28 and 130 M Q Myr -1 respectively. 
These rates are uncertain by a factor of a few, for several 
reasons, including the angle of the filament with respect 
to the plane of the sky. 

Although both accretion rates are uncertain to a factor 
of a few, it appears likely that the accretion rate onto 
the filament is larger than the accretion rate from the 
filament onto the central cluster. Formally, in order to 
have the latter rate be higher, the filament angle would 
need to satisfy 



28 M Q Myr- 1 130 M Myr 



tan(a) 



cos(a) 



(14) 



which is satisfied for a < 12°, i.e., a filament which is 
nearly parallel to the plane of the sky. Such a small incli- 
nation angle seems unlikely, but note that uncertainties 
in all observed quantities (mass, radius, velocity gradient, 
and infall velocity) could allow for a slightly larger range 
of angles. One must always use caution in interpreting 
velocity measures from different molecular tracers. The 
gas velocity could be some (non-constant) function of 
density, and using different molecular tracers for each ac- 
cretion rate measure could in principle reflect that, rather 
than overall flow differences in the two directions. The 
only way to circumvent this potential issue entirely would 
be to use observations of a single self-absorbed molecule 
and its rarer isotopologue for both measures. We do not 
have such data available at high enough signal to noise, 
however, the consistency of the HNC self-absorption dip 
with the N2H" 1 " and H 13 CO + peak emission suggests all 
are tracing similar volumes of gas, and that this should 
not be a major problem. 

If the tentative difference in accretion rates is true, this 
could imply that in the future, the accretion rate onto the 
central cluster will also increase (given the larger amount 
of mass available in the filament), and thus lead to more 
vigorous star formation in the cluster centre. It is also 
possible that the increase in mass in the filament could 
promote local instabilities in the filament, leading to core 
formation, and eventually to a series of star formation 
sites along the filament. Dense cores are, in fact, already 
observed in the AzTEC/ASTE data along the filament, 
as well as in 8 pm absorption. A small number of class 
and I YSOs are also seen, which may be an indication 
of past instabilities in the filament. The small cluster of 
class II YSOs just past the southern end of the southern 
filament suggests a previous episode of star formation, 
perhaps also associated with the filament. 

At the current estimated accretion rate of 
28 Mq Myr -1 , the cluster's membership would 
roughly double in 1 Myr with the addition of about 60 
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YSOs, assuming a mean mass of 0.5 Mq. An increasing 
accretion rate would lead to a larger number of YSOs. 
From these estimates, we can see that our estimated 
accretion rate from the filament onto the cluster is at 
a comparable level to previous observational estimates 
of the rate of mass supply required to form the current 
generation of YSOs in the cluster centre. We can 
therefore conclude that filamentary accretion is an 
important, if not the dominant, factor in governing the 
formation and evolution of stars in the central cluster. 

Throughout this analysis, we ha ve adopted a distanc e 
to Serpens South of 260 pc (e.g., iStraizvs et all 12003). 
for consistency with previous work, although recent 
VLBI observatio ns suggest that a distance of ~ 415 pc 
(|Dzib et all I2010D applies to at least the main Serpens 
cluster. If this distance is also appropriate for Serpens 
South, then the quantities we estimate in our analysis 
here become a factor of several larger. The mass in the 
southern filament would increase by a factor of 2.5, the 
velocity gradient would decrease by a factor of 1.6, and all 
other quantities of interest (lengths, accretion rates, mass 
per unit length) increase by a factor of 1.6; the accretion 
rates become 45 M Q Myr -1 parallel to the filament (onto 
the cluster) and 208 M Q Myr -1 perpendicular to the fila- 
ment. None of these changes effect our conclusions qual- 
itatively, and in fact, they tend to strengthen our finding 
that filamentary accretion is important in Serpens South. 

Filament-cluster systems such as Serpens South 
present a unique opportunity to gain a deeper under- 
standing of cluster formation. With more sensitive ob- 
servations, it should be possible to measure the spatial 
variation in the infall velocity across the filament, which 
could be used to model the full 3D velocity profile of ma- 
terial in the filament. This in turn can provide a straight- 
forward test of the relative importance of magnetic fields, 
gravity, and t urbulence in the ong oing cluster accretion, 
as discussed in lBalsara et all (|2001f ) and Heitsch (in prep). 

8. CONCLUSION 

Using a combination of dust column density infor- 
mation derived from AzTEC/ASTE continuum maps 
(Gutermuth et al, in prep) and dense gas kinematics de- 
rived from Mopra spectral maps, we study the flow of gas 
along the southern filament of the Serpens South embed- 
ded cluster. We find that the filament appears typical in 
terms of column density properties (e.g., a supercritical 
mass per unit length). The multiple molecular emission 
lines observed in our Mopra survey allow us to infer the 
possible presence of an accretion flow from the filament 
onto the central cluster (if the filament is in front of the 
cluster), and additionally we find that material continues 
to be accreted onto the filament. We estimate accretion 
rates of 28 M Q Myr -1 for accretion onto the central clus- 
ter, and a lower limit of 130 M Q Myr" 1 for infall onto 
the filament (assuming a distance to Serpens South of 
260 pc). Both of these are similar in magnitude to the 
estimated current rate of star formation in the central 
cluster. 

Filamentary accretion flows therefore appear to be an 
important mechanism for supplying the material neces- 
sary to form the YSOs in this cluster. As more embedded 
clusters and their associated filaments are identified in 
current surveys, it will be interesting to measure the as- 
sociated gas kinematics to determine whether other em- 



bedded clusters behave similarly to Serpens South. 
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Table 1 

Basic Line Properties 





Frequency 


HF components 


Ref. b 


Int. Range c 


Integrated Noise d 


Molecule 


(MHz) 


(if applicable) 




(km s — - 1 ) 


(K km s- 1 ) 


N 2 H+(l-0) 


93173.4669 


Fi,F = 2- 1,2- 1 


1 


-10.5-5-7.5 


4.38 


H 13 CN(l-0) 


86340.1840 


F = 2-1 


2 


-8.5-5-6.5 


3.87 


HCN(l-O) 


88631.8473 


F = 2-1 


3 


-9^7 


4.02 


HNC(l-O) 


90663.5560 


F = 2 -1 


1 


-2^2 


2.04 


H 13 CO+(1-0) 


86754.2884 


N/A 


5 e 


-2^2 


2.02 


HCO+(1-0) 


89188.5260 


N/A 


6 


-5^5 


3.23 


SiO(2-l) 


86846.9600 


N/A 


7 


-1.5-+1.5 


1.76 


NH 2 D(1,1) 


85926.2630 


N/A 


8 


-6.5-5-6.5 


3.66 



a In the case of a species with hyperfine splitting, the tra nsition use d for centering . 

b Reference for the adopted fre quency: 1. iPagani et all ()2009h : 2. iPearson et all (fl976l) . as recom me nded by JPL 
(http : / / spec . jpl . nasa . gov ) ; 3.fAhrens et~al (|2Q02f ), as recommend ed by CDMS (jMiiller et alll2005f) : 4. lBechtel et~al 
feooa): 5. [ Schmid-B ^rgketail (l20Ql, as recommended by CDMS; 6. lUlich fc Haai (fl976h. asrecommended by NIST 
( Lovas et a l 2009) 7. Mans on et all (jl977f ). as recommended by CDMS; 8. Main frequency listed for no hyperfine split- 
ting from Bes ter et all |l983), as recommended by NIST (our sensitivity is too low to pick up the satellite components). 
c The range in velocities, relative to the line centroid, used when calculating the integrated intensity. 
d The median formal 1-sigma error on the integrated intensity in brightness temperature units, i.e., the typical spectral 
rms multipled by the square root of the number of spectral channels integrated over. 

° Note that hyperfine splitting was detected by Schmid-Burgk04, but with the components separated by only 
0.133 km s _1 , this would not be visible in our Mopra data. 
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Table 2 

Hill5 Model Parameters Fit to HNC 



Fit a 


V , b 
' in fall 


8V h 


v lsr b 


r b 


n~< b 
x peak 




(km s - 1 ) 


(km s - 1 ) 


(km s - 1 ) 




(K) 


Fit-1 


0.25±0.02 


0.45 ±0.01 


0.05±0.01 


2.7±0.1 


6.33±0.06 


Fit-2 


0.54±0.01 


0.36 ±0.01 


0.06±0.01 


1.6±0.1 


6.6±0.1 



a In our 1000 randomized fits, two clusters of fit parameters were identified, as listed here. The original HNC fitted 
Hill5 parameters are all identical to the mean values listed for Fit-1. 

b Mean and standard deviations of the infall velocity, velocity dispersion (Gaussian sigma units), centroid velocity, 
optical depth, and peak excitation temperature, for the subset of the 1000 randomized trials. 
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Figure 1. The large-scale context of our Mopra observations of Serpens South. The greyscale image and contours show an AzTEC/ASTE 
map of the millimetre dust continuum emission, while the red filled circles, orange open circles, and yellow diamonds show the YSO 
population found in Spitzer data. The black circle at the bottom right corner shows the AzTEC/ASTE beamsize (28" FWHM). The blue 
contour shows the region mapped with Mopra, covering all of the main filamentary structure. 
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Figure 2. The observed emission lines near the cluster centre (-18 ,l 30 m 01 s , -2°03'12"). The data have been summed over 4x4 pixels 
(an area of 60 by 60") to increase the signal to noise level. The red vertical lines indicate the expected line centres based on a fit to the 
N2H" 1 " spectrum. In the case of (discernable) hyperfine components, all are shown, with the relative heights scaling with the optically thin 
line ratios. Note that SiO emission was not detected with certainty anywhere in the map. 
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Figure 3. The integrated intensity of the three brightest molecular tracers in our survey not affected by self-absorption (N2H" 1 ", H 13 CO + , 
NH2D), as well as SiO, which should not show self-absorption, but has no clearly detected emission. The greyscale images show the 
integrated intensities, with white at a value of -0.1 K km s _1 and black at a value of 90% of the maximum integrated intensities measured 
in each molecular species (23.5, 3.8, 3.2, and 1.2 K km s — 1 respectively). The grey contours in all panels are at values of 0.5, 2, 3.5, 7, 
and 15 K km s from darkest to lightest. The dashed line indicates the full region mapped by Mopra, and the blue contours show the 
dust emission from AzTEC/ASTE, at levels of 0.25, 0.5, 1, and 2 Jy bin -1 from light to dark. The Mopra and AzTEC/ASTE beams are 
shown on the bottom right of each panel in dark and light grey respectively. 
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Figure 4. The centroid velocity measured for N2H+ across Serpens South. The contours correspond to the dust emission levels shown 
in Figure [l] and [3] and the circles denote the same YSOs as shown in Figure [l] (dark grey corresponds to class 0, medium grey to class I 
and light grey to class II sources). The peak ridge of the filament is shown by the dashed lines (red in the south and blue in the north for 
clarity). The overall velocity gradient measured in the southern filament is indicated by the arrow at the bottom left, while the circle at 
the bottom right indicates the Mopra beamsize. 
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Figure 5. The toy cylinder model we consider for our accretion rate analysis. The cylinder has length L, radius r, mass M, and velocities 
of Vji along the long axis and Vj_ in the radial direction. It is inclined an angle of a (not shown) relative to the plane of the sky (with the 
top end further from the observer). 
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Figure 6. Left: the N2H+ (green) and H 13 CO + (yellow), and HNC (blue) average spectra for the southern filament. (For N2H+, only 
the isolated hyperfine component, (Fi,F=0, 1-1,2), is shown.) The HNC spectrum shows clear signs of self-absorption, while the N2H+ 
and H 13 CO" 1 " spectra are well described by a single Gaussian. Right: The filament-averaged HNC profile, along with the expected relative 
positions and strengths of the hyperfine emission lines of HNC. The hyperfine splitting is on too small a scale to be responsible for the 
observed self-absorption profile. The red line shows the best fit Hill5 infall model profile (dcVrics & Myers 2005|). See text for details. 
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Figure 7. A comparison of the self-absorption profiles for HCO + , and the three hyperfine components of HCN. The HNC emission is 
shown in dark blue, while H 13 CO + is in orange, the isolated hyperfine component of N2H+ (Fi ,F=0,1-1,2) is in green, HCO + is shown in 
light blue, and each component of HCN / H 13 CN is shown in light / dark pink. In the case of lines with hyperfine components, the velocity 
centroid has been shifted to the appropriate common centre. Note that no scaling has been applied to the brightness temperatures shown. 
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Figure 8. The radial column density profile of the southern filament in Serpens South. The solid black line shows the mean profile across 
the filament, while the grey shaded area indicates the standard deviation. The Az TEC/ASTE beamsize is sh own by the red dotted line. 
The light blue dash-dotted line shows the best fit profile, following the formulation in Arzoumanian ct al (2011), while the dark blue dashed 
line indicates the best fit profile when the power law index is that corresponding to an isothermal cylinder. 



